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Abstract 

The blowdown of high pressure gas in a pressure vessel produces rapid adiabatic cooling of 
the gas remaining in the vessel. The gas near the wall is warmed by conduction from the wall, 
producing radial temperature and density gradients that affect the flow, the mass efflux rate and 
the thermodynamic states of both the outflowing and the contained gas. The resulting buoyancy- 
driven flow circulates gas through the vessel and reduces, but does not eliminate, these gradients. 
The purpose of this note is to estimate when blowdown cooling is rapid enough that the gas 
in the pressure vessel is neither isothermal nor isopycnic, though it remains isobaric. I define a 
dimensionless number, the buoyancy circulation number BC, that parametrizes these effects. 
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I. INTRODUCTION 



The process by which a pressure vessel containing high pressure gas is vented to the 
outside through a narrow (compared to its diameter) orifice is called blowdown[l|. The 
rapid reduction in pressure and density within the vessel produces adiabatic cooling of the 
gas. If cooling is sufficient, gases of high condensation temperature may condense as liquid 
droplets, leading to a complex multi-phase flow. These droplets tend to fall to the bottom 
of the vessel, where contact with the warm wall will evaporate them, but if small enough 
may instead be entrained in the flowing gas. The blowdown rate will then depend on the 
position of the venting orifice and the nucleation process that controls the droplet size. 

In general, the gas inside the vessel may be regarded as isobaric because the blowdown 
time (for an orifice narrow compared to the vessel size) is much longer than the sound 
transit time across the vessel. Near the wall the gas is warmed by conduction from the wall. 
Typically, the vessel itself has a much greater mass and heat content than the contained 
gas, and may be regarded as a constant temperature heat bath. The thermal diffusion time 
across the gas is usually much longer than the blowdown time, so the temperature variation 
within the gas may be significant. 

Under isobaric conditions the temperature gradient implies a density gradient. In the 
presence of a gravitational field this density gradient drives a circulatory flow, with gas 
acquiring heat in a boundary layer at the wall and gradually mixing with the bulk of the 
gas. There is some similarity to Ekman pumping, but in the blowdown problem there is an 
additional imposed time scale, the blowdown time, that has no analogue in Ekman pumping, 
while Ekman pumping has a characteristic rotational velocity that has no analogue in the 
blowdown problem. 



II. TIME SCALES 

Here we estimate the relevant time scales and other parameters. The circulation time 
t C i rc is the characteristic time for gas to flow through the boundary layer in a (nominally 
spherical) vessel of radius r at speed v: 
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The boundary layer thickness is given by 



-Dtdrc, (2) 



where D is the thermal diffusivity. The circulatory speed v in a gravitational field g (or 
equivalent acceleration) may be estimated, approximating the flow as freely accelerated 
(approximately valid if the Prandtl number is of order unity as it is for dilute gases): 



where At is an Atwood's number: 



y/^At, (3) 



At « (4) 

tcooi 



where t coo i is the adiabatic cooling time. 

For unchoked flow of an ideal gas with adiabatic exponent 7 and sound speed c s from a 
volume V through an orifice of area A the adiabatic cooling time 

tco ° l = Ac a {n-iy (5) 

If flow through the boundary layer is fast compared to this time (t C i rc <C t coo i), as will usually 
be the case, At will be small because the gas in the interior of the vessel will only cool by 
a small fraction of its (absolute) temperature in the time required for gas in the boundary 
layer to flow through that layer. 

Combining the previous results we estimate the circulation time: 

tcirc ~ ( ] • (6) 

V 9 J 

The time required for the entire volume of contained gas to flow through the boundary layer 
is 

„ -7/6^/6 

f ~ f _ ~ cooi (7\ 

Lmix ~ t cjrc ~ /6 /2 . yi ) 
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If tmix > tcooi then the temperature gradients are large and it is not valid to treat the gas as 
isothermal or isopycnic (constant density). 



III. THE BLOWDOWN CIRCULATION NUMBER 

To facilitate dimensional insight, we rewrite the preceding equation: 

^mix ~ l cool l chaV \°) 

where the characteristic time, a quantity defined only by the properties of the vessel, gas, 
and gravity, is 

r 7/5 

tchar = gl/5£)3/5- ( 9 ) 

This permits defining a (possibly novel) dimensionless parameter, the blowdown circulation 
number: 

BC = r?/6 KIR « ^ « W 5/6 . (10) 

g^D 1 / 2 ^ tcooi \t CO oiJ 

If BC ^> 1 the gas is neither isothermal nor isopycnic and its circulation and density and 
temperature distributions must be considered in order to calculate the blowdown correctly, 
even if droplets do not condense. The escaping gas is drawn from near the surface of 
the pressure vessel, where it is warmed by thermal conduction from the wall, while the 
temperature in the central regions of the vessel drops nearly adiabatically. Because t mix is 
much greater than the actual blowdown time t coo i(j — 1), the gas issuing from the orifice late 
in the blowdown process will be much colder than its initial temperature. A quantitative 
calculation would, in general, require a three-dimensional numerical solution of the Navier- 
Stokes equations with heat flow. In the special case of an orifice at the top or bottom of the 
pressure vessel the flow would be two-dimensional with axial symmetry. 
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IV. PARAMETERS OF CIRCULATORY FLOW 



Explicitly, the flow parameters v, 5 and boundary layer Reynolds number Re are (for 
Prandtl number of order unity): 

l/3 r 2/3 
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For Re > 10 3 -10 4 the boundary layer flow will be turbulent and D should be replaced by 
a turbulent diffusion coefficient. Its value may be such as to reduce Re to the threshold of 
turbulent breakdown. Such large values of Re are only found for very large pressure vessels. 

As a numerical example, consider a spherical pressure vessel with r = 1 m with orifice 
A = 10~ 3 m 2 (10 cm 2 ) containing air at a pressure of 10 7 Pa (100 bar) and 20 °C. Then 
D = 2.12 x 10~ 7 m 2 /sec (0.00212 cm 2 /sec), t cooi = 31 sec, v = 0.68 m/sec (68 cm/sec), 
5 = 5.6 x 10~ 4 m (0.056 cm) and BC = 85. The time for gas to mix throughout the 
vessel t mix = BCt coo i = 2600 sec, much longer than the venting time V/ (Ac s ) = 12 sec. The 
temperature and density distributions within the vessel will be strongly inhomogeneous, even 
though the pressure remains uniform to high accuracy everywhere except in the immediate 
vicinity of the orifice, and, in general, three-dimensional computational fluid dynamics would 
be required for a quantitative calculation of the efflux rate and thermodynamic parameters. 



[1] This paper concerns pressure vessels that initially contain only gas. Blowdown of vessels that 
initially contain liquid or solid that flashes into vapor upon pressure release involves different 
phenomena. 
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